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ABSTRACT 

This paper is devoted to a study of possible scaling laws, and their logarithmic cor- 
rections, occurring in deep inelastic electropion production. Both the exclusive and semi- 
exclusive processes are considered. Scaling laws, originally motivated from PCAC and 
current algebra considerations are examined, first in the framework of the parton model 
and QCD peturbation theory and then from the more formal perspective of the operator 
product expansion and asymptotic freedom, (as expressed through the renormalization 
group). We emphasize that these processes allow scaling to be probed for the full am- 
plitude rather than just its absorbtive part (as is the case in the conventional structure 
functions) . Because of this it is not possible to give a formal derivation of scaling for deep 
inelastic electropion production processes even if one believes that they are unambiguously 
sensitive to the light cone behavior of the operator product. The origin of this is shown to 
be related to its behavior near a; ~ 0. Investigations, both theoretical and experimental, 
of these processes is therefore strongly encouraged. 
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1. INTRODUCTION 



1.1. Motivation 

The well-known logarithmic scaling violations in the structure functions of nucleons 
predicted by asymptotic freedom played a crucial role in establishing QCD as the accepted 
theory of strong interactions. These predictions, based as they are on the renormalization 
group and the operator product expansion of two electromagnetic currents near the light 
cone, are on a strong theoretical footing. This is in contrast to the situation in many other 
hadronic processes. For example, a theoretical analysis of the Drell-Yan process, or of jet 
production, requires some input of unknown soft infrared contributions . In spite of this 
ignorance, theorists concerned with these processes have gone ahead with recipes for their 
calculation, which often have been very successful. These recipes involve a careful mix 
of several ingredients, some of which such as asymptotic freedom and perturbation the- 
ory are well-understood, while others, such as the ingredient describing non-perturbative 
hadronization effects are chosen with an eye for the acceptability of the finished product 
rather than for their theoretical basis. Nevertheless this pragmatic approach is now taken 
to be sufficiently reliable that it is used to estimate backgrounds in experiments searching 
for unusual or new phenomena: for example calculations of jet production are used in 
searches for potential Higgs candidates. 

It is therefore important for theorists to appreciate the lack of knowledge which is a 
necessary input into such calculations involving perturbative QCD so that experimental 
results can be used to illuminate the approximations and assumptions used in the recipe. 
It is in this spirit that we consider here another hadronic process which is similar in 
many respects to the classic deep inelastic structure function process which underpins 
perturbative QCD. This is the process of deep inelastic exclusive pion production from 
nucleons which each of us worked on many years ago using the ideas of current algebra 
and PCAC. Our original analyses were clearly not rigorous but they did lead to predictions 
of scaling which agree with experiment albeit at relatively low energies. We therefore re- 
analyse this process from a modern viewpoint in this paper: we derive the scaling laws 
and calculate the expected logarithmic deviations from scaling. 

It seems to us that it is important to stimulate more interest in this and similar pro- 
cesses at this time because they are amenable to experiment in the not-too-distant future 
at HERA, CEBAF and a possible upgrade of SLAC. Much of what we say is applicable also 
to other processes which are described theoretically by an amplitude which is proportional 
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to the product of two currents. This should therefore include the deep inelastic electropro- 
duction of p's, K^s and most interestingly real photons: the latter process corresponding 
to a direct measurement on the non-forward Compton amplitude, albeit with one virtual 
and one real photon. We will show that if the predictions of scaling behaviour are indeed 
verified by experiment, important implications follow about the analytic structure of the 
amplitudes for these processes. 

In section 1.2 we define the process we are interested in and the corresponding kine- 
matic region. In section 2.1 we outline the derivation of the scaling law as was done 
originally by two of us over "twenty" years ago |^] ^ using Current Algebra. The result- 
ing predictions were, in fact, verified by some rather coarse data taken about that time 
0] P]. In section 2.2 we present a rederivation of the scaling law using the language of 
the parton model as done in the thesis of one of us in the early 1980's 0. In section 2.3 
we calculate logarithmic corrections to the amplitude using the diagrammatic approach 
pioneered in the seminal paper by Altarelli and Parisi 0. It should be noted that the 
end result of this approach is the integral for the amplitude M, which is analogous to the 
Altarelli-Parisi evolution equation. This equation can then be used to predict the moments 
of the amplitude, which in principle can be measured. The validity of the procedure is sub- 
ject to some reservations which we will discuss fully. In section 3 we approach the problem 
from the point of view of Operator Product Expansion. We explain why the prediction for 
the moment of the amplitude is sensitive to the analytic properties of the amplitude near 
x=0. This implies that an experimental study of deep inelastic pion production from which 
these moments can be determined may well yield information on the low x dependence of 
the amplitude. In section 4 we discuss the connection of our results to experimental data 
and suggest future experiments. 

1.2. Kinematics & Definitions 

The amplitude that we are going to study is defined as follows: 



= {p'7v\J^\p) 



(1.1) 




(1.2) 



(1.3) 
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Here J^(x) is the electromagnetic current, A^{x) the axial current, m.,^ the pion mass, 
/tt its decay coupling constant and 4>Tr{x) its field. In going from (1.2) to (1.3), the standard 
PC AC identification has been used: 

d^A^ix) = UmlcPAx). (1.4) 

The kinematics are illustrated in fig. 1: p is the 4- momentum of the struck target, p' its 
final momentum and q that of the virtual photon delivered by the scattered electron; q' 
will be used for the pion 4-momentum. 

The relationship to, and generalization from the amplitude probed by measuring the 
conventional structure functions is clear. In that case one is probing only the imaginary 
part of the forward Compton amplitude whereas in electropion production one measures a 
full amplitude which, in general, is non- forward. Formally, the difference can be expressed 
as probing the difference between a time-ordered, or retarded product, as in ( p..3|) , versus 
a commutator as in the structure function case: 

W^, = j d^xe''i-^p\[J^{x), JMM- (1-5) 

The full forward Compton amplitude is given by 

J^, = J d''xe''i--{p\T[J^ix),MO)]\p) (1.6) 

so that W^i, = Im J'^^. These are represented by the diagrams in fig. 2. 

A further crucial difference between the two cases is, of course, that in (|1.3| ), the 
kinematics of real pion production dictates that, even in the deep inelastic limit when 
q^ is large, q'^ must remain fixed at m^; on the other hand, in ( p..5| ), the magnitude of 
the virtual mass of both currents is always large in the deep inelastic limit. This latter 
condition ensures that the light-cone is unambiguously being probed and so justifies the 
use of the light cone operator product expansion. In spite of the fact that this is not clearly 
the case in pion production we shall argue below that a short distance operator product 
expansion may dominate the process when q^ is large. 

There is a subtlety in this procedure which is also present in the standard forward 
Compton amplitude case. The point is that this formalism leads to an expansion in powers 
of 1/x, where x = —q^/2p.q, and, in the physical region accessible to real experiments, |x| < 
1. Such an expansion therefore clearly does not converge. In the structure function case, 
this potential problem is finessed because, there, one is interested in only the imaginary 
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part, as in (L5), so an analytic continuation from the unphysical large (where the 
expansion presumably makes sense) to the physical region can be effected Q . Indeed this is 
why the results are expressed in the form of moments of the structure functions rather than 
the structure functions themselves. We shall discuss this in more detail below. However, 
it is clear from this, that in the pion production case is sensitive to a potentially 
interesting part of the formalism not readily accessible to the structure functions. Indeed 
it may well be that because of this "problem" pion electropro duct ion can cast interesting 
light on the general small x behavior as well as the general assumptions that underly the 
usual derivation. Before reviewing the old scaling arguments, however, let us recall the 
relationship between the measured cross-section and the matrix element M^: this is best 
done in terms of the tensor 

T^, = M^Mt. (1.7) 

The result is given by: 

(fa _ a E' iv^-qY^^ da 
dE'dQ'dn ~ 27r2g2 E 1 - e dQ' ^ ' ' 

where E{E') is the initial (final) energy of the electron in the Laboratory (LAB) system 
and u its energy loss: note that u = E — E' = p.q/M where M is the target nucleon mass. 
The polarization of the virtual photon is given by: 

6 = [1 - ^^'^^ ; tan^k]-^ (1.9) 
q Z 

where 9e is the electron scattering angle in the LAB. The quantity da/dQ represents an 
equivalent virtual photoproduction cross-section in the outgoing hadron center-of-mass 
(CM) system: 

da _ M^lg^lcM >t ^^-(T -T ^ 

dn 167r2Vf^2|g|^j^l2^ """^ ^^^^ 2'^ (1.10) 

- {qyu')eT,, + {-{2qyu^)e{l + e)}^T,,] 

Here W is the total CM energy so = s = {p + q)"^ = {j>' + q'Y ■ The z-axis is defined 
to be coincident with the direction of q whilst the electrons define the xy plane. Thus all 
of the (j) (azimuthal) dependence is contained in (T^a; — Tyy) ~ cos2(^ and Tn^z ~ cos(^. In 
what follows we shall limit ourselves to the case where the particle spins are unobserved. 

Finally, it is worth noting that in the deep inelastic limit {q^ —00), the square of 
the momentum transfer 

t^^^^{q'-qf = {p'-pf (Lll) 
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is constrained, in the physical region, to he between —2v and 

train ~ ~ T\ l-'---'-^) 

{\-x) 

Typicahy the hmit we wiU be considering keeps t and x fixed (and finite) with —oo. 
Thus X must not be too close to unity. Furthermore, the region of interest is predominantly 
forward scattering in the ttN CM system. In what follows it is convenient to write 

q^^En^ + A^ (1.13) 

where is a null-vector, i.e. = 0. Thus 

q'^ = 2E(n.A) + t (1.14) 

and 

-[2Ein.A)+t] 
i2[Ein.p) + A.p])- ^'-'^^ 

The scaling limit can then be realized by taking £^ — > oo with both x — (n.p)/(n.A) and 

t fixed. 



2. SCALING LAW 

2.1. Current Algebra 

Scaling laws for T^;^ can be derived using a current algebra approach augmented by 
some heuristic assumption about the light cone behavior of the commutator. This can be 
checked in perturbation theory and justified by the operator product expansion as sketched 
below. We begin by setting g'^ = in which case 

UM^^C^^q"" +E^ (2.1) 

where 

C^, = iJ d''xexp{iq-x)(j>'\eixo)[J^^ix),A,{0)]\p) (2.2) 

and 

E^ = J d''xexp{iq.x){p'\5{xo)[Mx),Ao{0)]\p). (2.3) 
Using the usual SU(2) x SU(2) current algebra 

S{xo)[Ai{x), J^m = Z6^3fe4(0)54(x) (2.4) 
6 



we immediately get that is independent of E; (it depends only on A and, from its 
usual parametrization, we get the well-known axial vector and induced pseudoscalar form 
factors of the nucleon). 

The scaling result we want to show is that C^^q''^ is also independent of E. We shall 
first sketch the derivation of this result based on the spacetime behavior of the current 
commutators. In what follows it is convenient to introduce standard light-cone coordinates 
for a four-vector as follows: 

a± = —{ao±az); a± = {ax,ay). (2.5) 

Then the scaling limit is equivalent to Q- ~ y/^i^ — > oo with q.^. ^ fixed. 

Causality allows us, at least naively, to replace 9{xo) in (2.2) by 9{x+) in which case 
the asymptotic behavior of C^i, is given by 

C^. ^ d^xexp{tq-x){p'\5ix+)[Ux),AiiO)]\) (2.6) 

Consequently, 

q'^C^, q-C^+ ^-J d^xexpiiq ■ x){p'\6{x+)[J^{x),AX{0)]\p). (2.7) 

The commutator in (2.7) can be expressed in the form 

iUx), A\{0)]6{x+) = A;,{x-)S\x^) (2.8) 

where A^{x-) is, in general, model-dependent. 

In QCD, canonical commutation relations lead to 

4(x_) = A;(0)5(x_)+S;(x_) (2.9) 

where i?^(x_) is an unknown (typically bilinear non-local) operator which is non-singular 
at X- ~ 0. We conclude then that, in the scaling limit, 

?'"C^. ~ j dx-exp{iq+X-){p'\B^{x-)\p). (2.10) 

This is the desired result since it shows that is independent of E and is only a function 
of X (through q+) and A. 

This result straightforwardly translates into the following scaling constraints on the 
components of T^i, occurring in the measured cross-section, (|1.8|): 

^{T^^ - Tyy) ^F2(x, t){kl - kl) ^2.11) 
T,, ^Fi(x,t)-F2(x,t)A2 

and Tj^z -F2{x,t)kj,/^z 
where A2 (as expressed in the LAB)^ — (|t + Muj) and the Fi are Lorentz scalars. 



7 



Parton Model 

We now turn to the treatment of this problem in terms of the quark-parton model. 
We assume that the two currents and A^, interact successively with the same quark 
while the rest act as spectators. Without QCD corrections the amplitude C^j^ is given 
by the sum of the two diagrams in fig. 3. This is analogous to the usual parton model 
treatment of the forward Compton scattering. 

To calculate these diagrams we assume that the struck quark carries a fraction r] of the 
momentum p of the hadron, in the infinite momentum frame. Immediately after absorbing 
the photon the virtual quark has momentum {r]p + q) , and hence it is highly off-shell. This 
is the basic reason for considering that the parton model is applicable here. The final 
quark has momentum {r]p + A) . Consider now the diagram in fig. 3(a). Its contribution 
is (tj are isospin matrices and Q is the generator corresponding to the electric charge): 

^ Tj A ^p'(r/P + A)75/t(?7^+ A)7mV^pM 
2^ 2[?7(n.j9) + (n.A)] 

To this has to be added the contribution from the crossed graph shown in fig. 3(b). A 
sum over the various types of quarks has been suppressed. V'p(^) represents the amplitude, 
or wave function, for finding a quark of a particular type carrying momentum k inside a 
nucleon moving with momentum p. The complete matrix element requires an integration 
over ?7, consistent with the requirement that {r}p + qY > i.e. rj > x. Schematically, the 
parton contribution is thus given by 

Ml = [ dr]iJp,{r]p + A)M^{r],p,A)i;p{r]p). (2.13) 

J X 

The matrix Al^ can be read off from (2.13) with an additional contribution coming from 
the crossed graph. This shows explicitly that, when E is large, M depends only on t and 

X. 

Note, incidentally, that for the parton model description of the conventional structure 
function, A = and only the imaginary part of M, i.e. its delta- function contribution, 
is required. In that case the integral reduces to xf{x), where /(ry) = il>p{r]p)i^p{r]p), 
the probability for finding a quark with fraction rj of the total momentum. It is worth 
remarking that in the analogous kinematic configuration here, where A f^i 0, the full Born 
amplitude reduces to 

M^^[Q,^]{n.p)p^. (2.14) 
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2.3. Leading Logarithmic Corrections 

In this section we sketch a computation of the leading logarithmic corrections to 
the parton model result. Rather than give a complete detailed description we present 
here only the salient features for the simpler and more limited kinematic situation where 
A 0. The point is that a subgraph which gives a logarithmic contribution ~ InA^ for 
A 7^ obviously is singular in the A ^ limit: hence the same diagrams (i.e. ladders and 
self-energy insertions to external legs) that give the leading logarithmic corrections in the 
usual deep inelastic scattering case (A = 0) also give the leading logarithmic corrections 
in the A case. 

For ease of presentation, we shall from here on use p to denote the momentum of the 
initial quark rather than of the initial hadron. The initial quark is taken to be off-shell by 
an amount comparable to the inverse of the confinement radius of the nucleon. Since the 
momentum transfer A is typically of the same order of magnitude, we make no distinction 
between \n{—q'^/p^) and ln(— qf^/A^). When, for example, we end up with a logarithmic 
integration of the form J dk'^/ {k^ + A^) we shall be free to take the A — > limit and write 
it in the form J^a dk'^/k'^. Although the final scattered quark is slightly off-shell, the A — > 
limit enables us to take the 7 matrices between on-shell spinors. It is important, of course, 
to take the A — > limit after having secured that the final integration is logarithmic. 

It is convenient to work in the light cone gauge where the propagator for a gluon of 
momentum k is 

G^Ak) ^ 1^ (2.15) 

with 

D,.{k) = g,.-^^^^^^^^. (2.16) 

In such a gauge the only diagrams (apart from self energy parts) giving leading logarithmic 
corrections are the ladder diagrams of fig. 4(a). It may be mentioned, in particular, that 
diagrams of the type of fig. 4(b) do not give leading logarithmic contributions. Incidentally, 
had we followed the operator product expansion approach these diagrams would correspond 
to contributions from gauge non-invariant operators. 

Turning now to the calculation of the ladder diagrams we show that a familiar picture 
emerges: there is strong ordering in the momentum flowing through the ladder and an 
evolution equation can be derived. We first examine the contribution from the lowest 
order diagram (fig. 5). Using a Sudakov parametrization for the quark momentum k 

k = ac + Pp + kj_ (2.17) 
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(2.18) 



(where s 



/x). We obtain: 



ri -agCp s f ,,9 , ,^ 



(2.19) 



7^(/:+ A) A(/fe+ A)7m /^7p75 



D 



[{k + En + A)^ + ie] (A;2 + ie) [{k + A)2 + ze] [{p - k)^ + ie] 



Initial and final spinors have been suppressed. Note that the amplitude is color singlet 
in the t-channel and that, in terms of Sudakov variables, 

k^ = a(3s-ki, {p - kf = -a{l - p)s - kl, {k + Af = aps - k^ - asx. (2.20) 

We first perform the a integration. In the region x < P < 1 there is one pole at a = 
—k'j_/{l — P)s due to {p — kY lying below the real axis whereas in the region Q < P < x 
there is one pole at a = k\/ Ps due to k"^ lying above it. In both cases we close the a contour 
so as to pick the contributions from those poles. In the regions /? < and /? > 1 all the 
poles with respect to a lie on one side of the real axis and can be avoided (we do not take 
{k + En + A)^ into account since it will be combined with the next element of the ladder). 
Observe that the leading logarithms come from the wide range of integration <k\<^ s 
where ji is an arbitrary renormalization scale. The parameter P is finite (typically of order 
x) whereas a ~ —k\/s is small. Hence the logarithms come, as expected, from the colinear 
configuration. 

After the a integration the denominator behaves like {k\Y, so we have to extract 
one k\ from the numerator if the final integration is to be logarithmic. Having done this 
we can pass to the coUinear configuration k = Pp, o; fa 0. As already remarked we shall 
also set A ^ and take the 7 matrices between on-shell spinors u{p). One might worry 
whether in the limit A ^ we lose logarithms multiplied by q.A/q^. There is no such 
danger since we parametrize everything from the start in terms of vectors n, p and A. Then 
q.A/q^ ~ E{n.A)/2E{n.A) ~ i. There will remain a factor 75/1/^7^ from the numerator 
which will combine with {k + En + A)^ from the denominator to form the Born term 
Mb{Pp, k'j_). Note finally that the logarithmic contribution comes from the region P > x, 
so that the propagating quark line remains highly ofi'-shell. 

The numerator in the integrand has the form (apart from the 75): 



Cp{p-k)a + c^{p-k)p 
c.{p — k) 



}■ 



(2.21) 
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After some algebra this can be reduced (in the A hmit) to a famihar form for the 
leading lowest order correction: 

dkl as ^ , 2/? 



M^{x,q')=2CF -gr^ d^{l-^ + ^-^}MB{x/f3,ii') (2.22) 

Iterating this an arbitrary number of times leads to an evolution equation: 

M{x,q^) = M{x,i?)+2CF / d(5^^M{x/[5,kl). (2.23) 

Up to now we have considered skeleton graphs only. When we dress the ladder with 
vertex and self-energy corrections further leading logarithmic contributions coming from 
the ultraviolet region are induced. These can be taken into account simply by replacing 
the "bare" coupling constant ctg in ( p.23| ) by the running coupling constant as{k\) ~ 
47r/(/31og/c^) where (3 = (11 — 271//3) (n/ being the number of flavors). In addition, 
the second term in ( |2.23| ) must be multiplied by the quark wavefunction renormalization 
constant Zp in order to cancel the infrared divergences from the soft gluon region and get 
a gauge invariant result. Zp in the light-cone gauge has been calculated in a number of 
places. The final result is 

M{x,q^)=M{x,ii'')+2CF ^^^^^ d(3P{(3)M{x/(3,kl) (2.24) 

Jf,2 4n 

where, 

P{(3) = - 5(1 - P) dxj^. (2.25) 

Equation ( |2.24| ) is somewhat difficult to handle from the phenomenological point of 
view. If for the moment we disregard the subtleties regarding the low x dependence of M 
which will be discussed in the following section, then we can disentangle M in ( 2.24 ) by 
taking moments. Defining 

Mniq^)= [ dxx^'Mix.q^) (2.26) 
Jo 



we get 

M (n'^\ Tin ^ 

where 



M^{q^) ^ (ini.Y'^ (2.27) 



rf^ = % f 1 + 4 y - , 1 (2.28) 

/5 V ^ J (n + 2)(n-F3)/ ^ ' 

3 — ^ 
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3. OPERATOR PRODUCT EXPANSION 

In this section we discuss an operator product expansion (OPE) analysis of this am- 
pUtude. As aheady intimated, there are some subtleties that prohibit a straightforward 
prediction for its asymptotic behavior. Before discussing this, however, it is worth re- 
viewing briefly the standard treatment of the conventional structure functions from this 
viewpoint. In that case the use of perturbation theory to determine the large behavior 
can be justified from the application of the OPE to the light cone expansion of J^v, (|1.6|). 

Explicitly, the time-ordered product of the two currents can be Taylor expanded 
around ~ in terms of a complete set of operators O^^^ "^": 

r[J(x)J(0)] ^ 5^c^(a:2)a;^^....a:^„O^;,--'^"(0). (3.1) 

[For ease of presentation, the currents are here taken to be scalar] . Equivalently, its Fourier 
transform is given by 

/ dW^-T[J(x)J(0)] ^Y.^mn{q^f-ff^O^^-^- (3.2) 



where 



Cr^^iq') ^ ( - I rfW^-^c^(x2). (3.3) 

On dimensional grounds the Cj^nio^) behave, up to logarithms, like {q^)~'^'' for large 
q^ where dc = 2dj — 4 — {do — n) is the dimension of C, dj that of J(x) and do that of 
the operator Ol^;^"^" . This is, of course, the origin of the observation that the asymptotic 
behavior is controlled by the operator having the lowest twist tq = do — n. Notice that 
these equations are all properties of the current operators (i.e., in the case of interest here, 
the pion and the virtual photon) and do not depend on the target state. For forward 
scattering we require the ground state target matrix elements 

{P\0'^'n-'''^\P) = AmnPf,^--Pf^^ + Bmn9 f,,P f^, ■ ■ ■ ■P^^r^ + ' (3.4) 

The Amm Bmn etc. are simply numbers characterizing the target. In the contraction 
of this with (^) it is clear that, in the Bjorken limit, terms involving the Amn dominate: 



one thereby obtains 

J{x,q^) = j <iW"'{p|r[J(l)J(0)]|p> 



E 



AmnCmni^q 



(3.5) 
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The large behavior of the CmniQ'^) can be determined from the renormahzation group 
using the asymptotic freedom property of QCD. Typically, for the leading twist operator, 
the Cjnn are dimensionless and behave like (In q^'^~°'^-^ where amn is determined by the 
anomalous dimensions of the Omn- Implications for the structure functions, which are the 
absorptive part of J , can be obtained using the standard analytic properties of J . This 
leads to the well-known result relating the moments of W to Cjnn{,Q^)'- 

M^{q^)= j dxx''-'^[vW{x,q^)] 

m 

The sum over m is, of course, finite and typically contains only a rather small number 
of terms. QCD therefore gives a specific prediction for the g^-dependence of each moment 
and it is this that has been successfully checked against experiment |^]. 

Now, suppose that experiments could be performed that directly measure the large- 
q^ behavior of the full amplitude J{x,q^). What is the QCD prediction for this? One 
immediately sees the difficulty: the expansion, ( |3.5|) , presumably only makes sense for 
I a: I > 1 and this is outside of the physical region. Indeed the analytic continuation to 



|x| < 1 ultimately leads to the moment equations, (|3.6| ) . Ideally, one would like to 
have a complementary expansion valid for |x| < 1; this would require knowledge of the 
analytic structure near x ~ which, unfortunately is not reliably determined by the 
RG. Naively, one could proceed with J just as one proceeded with the M^; i.e., simply 
take ^ oo in (|3.6| ) and pick out the dominant Cmn{q^) as determined by the smallest 
anomalous dimension. In the singlet case, for example, the conservation of the stress- 
energy tensor means that it has no anomalous dimension and so M2{q'^) asymptotically 
approaches a constant. This, in turn, means that the leading behavior of the Ti^2(?^,a:), 
the two conventional scalar amplitudes occurring in the decomposition of J^^^, is given by 

^,(,^.)«ZM^«£«!^f,^) (3.7) 



2x nq'^x \16-f3n 

Now let us examine the extension of this to the non-forward case. Eqns. ( p.l|) - ( |3.3| ) 
remain valid since they are properties of the currents and the expansion ( p.l| ) is supposed 
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to be in terms of a complete set of operators; (|3.4| ) however, clearly needs to be generalized. 
This can be straightforwardly accomplished by writing: 



fc=0 

+ Bjnjik{t)g^^^^p^^ ■ ■ ■ p^f^A.^^^-^ ■ ■ ■ A^^ + ■ ■ ■] 

Clearly ^^^^(0) = Amn and Bmnn{0) = Bmn- When contracting this with ( |3.2| ) we 
shall need the quantity: 

2A.Q t , , 

—i = \^-. 3.9 

It is, therefore, the wider set of coefficients Amnk{i) that dominate the asymptotic behavior: 
( P?5|) is thereby generalized to 

J{x,q ,t)^ . (3.10) 

n=0 m k=0 

The Cmnio^) are the coefficients appropriate to the axial current case of interest here, 
as expressed in ( p..3| ) and (|2.2| ) , and are the analogs of the Cmn{q^) of ( |375| ) . They, 



too, generally fall with q^ like powers of In(g^) determined by the appropriate anomalous 
dimension. PCAC ensures that, in M^, there is an operator with vanishing anomalous 
dimension so that it becomes a function of q^ and t only. In any case, the corrections to 
this will, as usual, be powers of ln(— g^). As already explained it is not possible, beyond 
this, to give the precise prediction for the large q^ behavior without summing the series. 



Finally, it should be noted that the result expressed in (|3.10|) is clearly not valid unless 
t <^ q^, which means that x must not be too close to 1. On the other hand, probing scaling 
and its violation should shed some light on the a; ~ region: if the predictions of this paper 
are experimentally verified it means that the relevant amplitudes are smooth in the region 
of small X where it would appear that the operator product expansion breaks down. 



4. COMPARISON WITH EXPERIMENT 

In this Section we begin by reviewing the connection of our results to the existing 
experimental data. At present, the only such data is for the inclusive reaction (and this 
was taken over 20 years ago at Cornell). It is possible, however, to extend the above 
arguments to this case provided the mass of the final hadronic "target" state {W) remains 
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(fa 

s- — ■ ^ F{x,t,W'^). (4.2) 



relatively small. The main difference is that the scaling function will now depend on W' 
in addition to x and t so that, instead of 

s'^^F{x,t) (4.1) 

which is the result implied by (|2.11j ) for the purely exclusive case, one now expects for the 
inclusive 

dtdW 

In the Cornell experiment, d'^a/dtdW'^ was measured at two different values of 
v^(2.66 and 3.146*6^) but at the same value of x. The data was averaged over 9 and 
(j). The scaling result, (^]^) , implies that the spectra, when plotted as a function of W , 
should be identical apart from a normalization factor (3.14/2.66)^ ~ 2.41. The data, as 
can be readily seen in fig. 6 are in remarkably good agreement with this prediction. 

These data are also presented in terms of the transverse momentum P± , the transverse 
momentum of the pion relative to the direction of the incoming virtual photon and of a 
variable x' which depends on the longitudinal momentum of the pion: 

x' = — r (4.3) 

Here Pn is the pion momentum along the direction of the virtual photon and Pmax is the 
maximum pion momentum. In fig. 7 

(4.4) 

is plotted as a function of P]_ at the two values of W mentioned previously and at two 
different values of x' . The straight lines are fits of the form Aexp{—BPj_). The similarity 
of the spectra suggests that the P]_ distribution (for fixed x) does not depend on g^, again 
in striking agreement with the scaling argument. 

It is clearly important for new experiments to be carried out at HERA on deep inelastic 
pion electroproduction at high energies check whether scaling continues to hold subject 
to the logarithmic violations which follow from ( |2 . 24| ) Furthermore the arguments of this 
paper do not only apply to deep inelastic pion electroproduction. A similar argument could 
be made for deep inelastic electroproduction of any particle which couples to a nucleon by 
means of a local current operator. This would include deep inelastic electroproduction of 
real photons, or of lepton pairs, or of p's, or K^s, or i/^'s, or T's, for example. It would 
be especially interesting to check whether the amplitude for the production of each of 
these particles scales in the same way, or whether processes which involve heavy quarks 
are different. 
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5. CONCLUSION 



In this paper we have shown how scahng laws for deep inelastic electropion production 
derived on rather general grounds from QCD-inspired current algebra, are manifested in 
QCD perturbation theory. Leading logarithmic corrections are calculated and an evolution 
equation for the amplitude derived. These are quite similar in character to the well- 
known ones occurring in the conventional Compton amplitude but have the advantage 
that the predictions for the full amplitude are, in this case, amenable to experiment. In 
the Compton case, only the imaginary parts (the conventional deep inelastic structure 
functions) are, in practice, measurable. However, for the full amplitude we show that, 
contrary to one's naive expectation, the usual deviations from scaling derived from an 
operator product expansion analysis do not lead to a well-defined prediction in the physical 
region. Thus, unlike the structure function case, the QCD perturbation theory result 
cannot be "rigorously" justified from asymptotic freedom. The reason for this can be 
traced back to the behavior of the amplitude near x ~ 0; the OPE leads to an expansion 
in l/x which cannot converge for a physical process. The conventional moment equations 
for the structure functions which exploit the known analytic properties of the amplitude 
are precisely designed to circumvent this difficulty. Thus, observation of the scaling laws 
and their violation for the full amplitude can potentially shed light on the small x behavior 
and help clarify just how far one can push results based on QCD perturbation theory. 

With renewed interest in such problems stimulated by recent HERA results and the 
potential of detailed data from CEBAF (albeit at relatively low energies) we feel that it is 
important to examine processes such as these that are natural extensions of the canonical 
structure functions. It should be stressed that these processes should also be viewed as 
yielding complementary data on the quark-gluon structure of the nucleon. In future work 
we intend to explore this aspect of the problem in more detail; meanwhile, the main thrust 
of this paper is motivated by the desire to rekindle interest in such problems. In fact it 
was partly stimulated by a query along these lines from the experimentalist Bogdan Povh 
of the University of Heidelberg; one of us (GBW) would like to thank him for his original 
inquiry. 
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6. FIGURES 




the external particles. 




conventional structure functions, ( |1.5|) . 
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crossed contributions. 
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Fig. 6 

The virtual photoproduction cross-section at two different values of W = ^/s but at 
the same value of u = 1/x; the data are taken from 0. According to eqs. (4.1) and 
(4.2) these should be identical except for a scale factor W^. 
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Fig. 7 

Transverse-momentum distribution of the produced pions for two regions of longitu- 
dinal pion momentum. The similarity of the data at different values of W are in 
agreement with the scaling argument. 
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